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Recent developments on higher dimensional gravity resulted in a number of interesting theoretical ideas such as 
the brane-world concept in which the standard model fields are confined to our four-dimensional world viewed as a 
hyperspace embedded in the higher-dimensional spacetime where gravity can propagate [U, 0, The simplest models 
in this context are proposed by Randall and SundrumQ. It was argued that the extra dimensions need not be compact 
| and in particular it was shown that it is possible to localize gravity on a 3-brane when there is one infinite extra 
. dimension [H, 0| ■ A striking consequence of the theory with large extra dimensions is that it can lower the fundamental 
\& ' gravity scale and allow the production of mini black holes in the universe. Such mini black holes are centered on the 
, brane and may have been created in the early universe due to density perturbations and phase transitions. Recently 
it was proposed that such mini black holes may also be produced in high energy collisions and could be probed in 
particle accelerator experiments in the very near future 0, [H, 0, 0|- 

A great deal of effort has been expanded for the precise determination of observational signatures of such mini 
black holes. One among them is the black hole quasinormal modes (QNM), which can disclose the evolution of the 
perturbation around black hole during a certain time interval and was argued carrying unique fingerprint of the black 
hole existence and is expected to be detected through gravitational wave observations in the near future (see reviews on 
this topic and references therein Q). Recently the QNMs of a brane-world black hole have been studied in [jjl Hoi. fll| . 
The gravitational perturbation of the mini black hole would be a characteristic sound and could tell us the existence 
of such black hole. Another chief possibility of observing this kind of mini black hole is the spectrum of Hawking 
radiation which is expected to be detected in particle accelerator experiments [ll|, EH EH EH EH EH EH El; El; H3| • 
Since such small black holes carry information of extra dimensions and persist different properties compared to those of 
ordinary 4-dimensional black holes, these two tools of detecting mini black holes can help to read the extra dimensions. 

Most of the available works have been done for the idealized case by studying the brane-localized modes in the QNM 
and Hawking radiation. One of the reason is that in the brane world scenario it was assumed that most standard 
matter fields are brane-localized. In the study of Hawking radiation, it was argued that the emission on the brane 
is dominated compared to that off the brane [2 1[. This argument was sup por ted numerically in the case of standard 
model field emission by the (4 + n)-dimensional non-rotating black holes [l7l . [Lsj and also in the higher-dimensional 
rotating black hole background . Some counterexamples to this argument have also been found . However in 
general it is very hard to obtain exact soltions of higher-dimensional Einstein equations so that the knowledge of 
the bulk is lacking, the thorough study on the field emission in the bulk is still a challenging task. Using a recently 
constructed exact black hole localized on a 3-brane in a world with two large extra dimensions [22| , Dai et al explored 
the Hawking decay channels with the influence of the brane tension [2(|. They found for the non-rotating black holes 
the dominate channels are still the brane-localized modes. 

In this paper we are going to study a braneworld black hole obtained in (23| . A general class of spherically symmetric 
and static solution to the field equation with a 5-dimensional cosmological constant can be derived by considering a 
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general line element of the type 



ds 2 = A(r)dt 2 - -—dr 2 - r 2 dfl 2 , (1) 
B{r) 



and relaxing the condition A(r) = B 1 (r) used in obtaining the usual 4-dimensional black holes. Casadio et al 
obtained two types of solutions by fixing either A(r) or B(r) and in one case with the choice A{r) = 1 — 2M/r, the 



metric reads [23 



1M 1 - — 
ds 2 = (l-_)(ft 2 d r 2 -r 2 dQ 2 . (2) 



This braneworld black hole is called CFM black hole. It would be noted that Eq.(2) was not derived from 5-dimensional 
solution, but rather was obtained as a solution to 4-dimensional equations constraining the possible form of the 5- 
dimensional metric. The extension of the asymptotically flat static spherically metric on the brane into the bulk has 
been discussed in [24|. The usual 4-dimensional Schwarzschild black hole is recovered with 7 = 3. The corresponding 
Hawking temperature is given by[23j Tbh = yl 3(7 — 3)/2/ (8irM). Thus in comparison with the Schwarzschild 
black hole, the braneworld black hole will either be hotter or colder depending upon the sign of 7 — 3. In this work 
we are restricted to the case when 7 < 11/3 as in |10| to ensure the Hawking temperature to be physical. The QNM 
of this kind of braneworld black hole was studied in [1J]| ■ We will concentrate on the absorption problem of this black 
hole. Immersed in external radiation fields, the black hole scatters the radiation and may either absorb or amplify 
it. The process is described by dynamical perturbation equations in the black hole background. We will study the 
scalar and axial gravitational perturbations on the background of this CFM black hole and investigate the absorption 
problem. We will extract the information on the extra dimensional influence on the absorption of the braneworld 
black hole by comparing to that of the usual four-dimensional black hole. 
For future convenience, we adopt the notation in 25], 

ds 2 = e 2v dt 2 - e 2 ^ 2 - e 2 ^dr 2 ~ e 2 ^d9 2 . (3) 

Comparing with (2), we have 

e 2 » = A(r), e- 2f * 2 = B(r), e 2 ^ = r 2 , e 2 ^ = r 2 sin 2 9. (4) 

In the vicinity of the event horizon rjf, the asymptotic expressions of the metric components A(r) and B(r) can be 
written in the first order 

A(r) ~ A 1 {r/r H - 1), 

B(r)~Bi(r/r ff -l), (r^r H + 0), (5) 

where A\ — 1 and B\ = 4—7 arc dimensionless, positive real numbers. Demanding that the spacetime is asymptotically 
flat, we have, 

A(r) -> 1, 

B{r) -> 1, (r->oo). (6) 

Our paper is organized as follows: in §11 we will go over scalar perturbation and investigate its energy absorption; 
in §111 we will derive the wave equation of the gravitational perturbation and study its energy absorption spectra. 
We will sumarize our results in the last section. 



II. SCALAR PERTURBATION AND ITS ENERGY ABSORPTION 



For simplicity we first consider the braneworld black hole immersed in the massless scalar field confined on the brane. 
We will go over the scalar perturbation and study the energy absorption problem in the general four dimensional 
spherical metric, Eq. ((T|) and then apply the general results to the CFM black hole. 
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A. Scalar Perturbation Equations 

The scalar perturbation is governed by the Klein-Gorden equation. For the massless scalar field, we have 

= i~9)~ l/ % [{-g) 1/2 g^d^] = 0. (7) 
Using the decomposition of the scalar field 

oo 

* =J2 R i( r ) p i ( cos ) e ~ tut > ( 8 ) 

1=0 

we have the Schrodinger like radial wave equations, 

-^-l+[^-V s , l (r)](rR l ) = 0. (9) 
r* is the tortoise coordinate defined as dr* = dr j\J AB. The effective potential reads 

V Sil =A l -^A + ^(A'B + AB'). (10) 

B. Boundary Conditions 

There are freedoms in choosing boundary conditions depending on physical pictures of different problems. We will 
consider the plane wave scattering in our work. 

1. Boundary Condition at Infinity 
In the remote region, by virtue of the scattering problem, the scalar field at infinity is in the form 



gicjr cos _|_ f iy) giLor 



e~ luJt . (11) 



Referring to the scheme of partial wave method in quantum mechanics, the plane wave can be expanded into spherical 
waves as, 

oo 

e iurcosd = ^(2Z+l)i'ji(wr)Pi(cos0). (12) 



1=0 

When r tends to infinity, 



27-1-1 
2ior 

1=0 

The scattered wave can also be expanded in the form 

Me- = f> + l)/f^P ; (cos0), (14) 
r 2tor 

1=0 

where /* are complex constants, and the superscript "s" indicates the case of scalar field. Inserting this equation 
together with Eq. (fl~3"|) into Eq. (fTTj) , we have the boundary condition at infinity 

°° 97-1-1 

$ „ e -^ty^(_l)i+i_^[ e -^'-. _ (-iysf e ^«m(cos0), (15) 
^— ' 2ur 

1=0 

where Sf = 1 + /f, and the radius r on the exponentials are replaced by tortoise coordinate for future convenience. 
Note that at infinity c?r* ~ dr so that such replacement on the exponentials makes no difference. In accordance, we 
have the boundary condition of Ri at infinity 



91 4- 1 

Rl ~ ^—[ e -^ r * - (-lYS!e iur *}. (16) 

2u;r 
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2. Boundary Condition at Horizon 

Now we examine the boundary condition in the near-horizon region. We can express the solution of Eq. ([9]) in 
power series [27j 

oo 

Ri = (r/r H - 1)? di, n (r/r H ~ !)"• (17) 

n=0 

In the vicinity of the horizon, we have the asymptotic behavior 



VAB ~ y/AiBi(r/r H - 1), (18) 

as well as 

V s ,i(r)~0{r-r H ), (r^r H + 0). (19) 
Inserting the expansions Eq. (I18I19|) into the radial equation Eq. ([9]), we have 

A 1 B 1 \p(fi-l)+p]+r%u i = 0, (20) 

or, 

Since nothing can escape from the black hole, only the minus sign is acceptable. Thus we have the boundary condition 
of Ri(r) near the horizon [27| 



Rl(r) ~ 9i(r/r H - 1)-*-^/^^, (22) 

where gi are the zero-th coefficient of the expansion of Rf. gi = d^o, which is dimensionless. The boundary condition 
for the total perturbational scalar field $ is given by 

oo 

$ rsj e ~ luJt ^gi{r/r H ~ ir MrH/v/7 ^Pi(cos6l). (23) 

1=0 

C. The Energy Flux and the Absorption Cross Section 

The energy flux is derived from the energy-momentum tensor of the scalar field. For massless scalar field, it is 
generally expressed as [28) 

T av = - l -g^ a ^ a . (24) 

The energy-momentum flow of the spacetime is defined by[26j 

P»=Tj??(t), (25) 

where vector £(t) = d/dt is the time- like Killing vector. We now calculate the energy falling into the black hole per 
unit time using Gauss' theorem. In the theorem, for any appropriate four vector and space time sector AT, there 
holds 



A»flt= / A^Os, (26) 
N ' JdN 

where r?4 and $^3 are the volume elements of the spacetime sector Af, and its boundary dAf. n M is the normal vector 
of the boundary dAf which satisfies n A = f^, where n is the 1-form physically equivalent to vector n M . Here we 
may choose the space time sector Af to be the spacial area within two spheres of radius T\ and r 2 , and with the time 
interval between t\ and t2, whose boundaries are dAf. To find out the energy falling into the hole, we focus on the 
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hypersurface of r = r\, whose normal vector is = (0, vM7B,0,0), and we let r% — > rjy. Thus the energy falling 
into the hole within the time interval between t\ and £2 is [25:], 



E {ahs) = / di / P»n fi r 2 H sin9d6d(p= dt T ^r 2 H sin 6 d6d<p. 

Jt\ J Aix J t\ J At: 

Applying the explicit forms of T 0fl and n M , as well as the boundary condition Eq. (p£5|) . we have, at the horizon 

00 

9i Pi (cos ( 



(27) 



z=o 



Inserting this into Eq. (|27p , we have the total energy falling into the horizon per unit time: 

00 



Jfp(abs) 00 |„ 12 

$ (abs) = dh ' a_. .2_2 Iff I 



Here, we have used the orthogonality of Legendre polynomials, 



1=0 



21 + 1 



Pl(x)Pl,(x)dx . )/+| . 



(28) 



(29) 



(30) 



On the other hand, following the similar procedure as above, we can derive the expression of the energy flux of the 
incoming wave at infinity, which is simply 

Therefore for the massless scalar wave the total energy absorption by the black hole is given by 

§(abs) |„.|2 



r abs 



(inc) 



1=0 



21 + 1 



For angular index /, the partial absorption cross section is expressed as 

13/ 1 2 



af, s = 4irr 2 



H 



21 + 1' 



After we get gi we can obtain the final result of the absorption spectrum. 

For the radial equation Eq. ©, the Wronskians for any two solutions i£j and rS 2 ^ is 



W[R^,R^] 



r 2 VAB 



(32) 



(33) 



(34) 



with the constant K12 to be determined by the explicit form of the two solutions. Using the asymptotic solution of 
Ri at infinity Eq. (fT6|) . we have 



W[RIR^ 



i(2l + lf 



1-13 



s\2\ 



2Lur 2 VAB 

and employing the asymptotic solution at the horizon Eq. (|22p . we get 



W[RlR l ]=-^L\g l \ 2 



r 2 \JAB 



Equating the above two equations, we find 



1.9, 



(21 + l) 2 
Aijj 2 r 2 H 



(1-131 



s\2\ 



Thus with Eq. p3[) we can express the energy absorption in the form 

7T 



abs 

a, I ~ ~2 



(21 + 1)T 



1 1 



where Tf is defined in terms of S* ; s as 



r? = 1 - I3T- 



(35) 
(36) 
(37) 

(38) 
(39) 
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Figure 1: The effective potential of the scalar perturbation with the change of 7 and I. 

D. Numerical Study for the Scalar Perturbation 

Now we proceed to calculate numerically the coefficients gi. Instead of solving Eq. ([9]) directly for Ri, we choose 
the particular solution ipi (r) , normalized as 

y l {r)=gi 1 R l {r). (40) 
In the asymptotically fiat region, it resolves into outgoing and ingoing waves 

W (r) = //- ) ^(+)+// +) ^(-). ( 41 ) 

where H E^j 

n+ I 

\ r J 

n=0 

satisfying = ^it-) an d the coefficients rj^ being unity for simplicity. The combination coefficients / ; are 

called jost functions. Comparing this to Eq. (fTo) . we have 

(-l)^ 1 / + 1/2 

9i = rrj - ' (43) 

w ft ] 
A+) 

St = (-l/ +1 4zy. (44) 

Thus if the jost functions are known, we can get the solution to our problem. Applying Eq. (|34|) to the two particular 
solutions at infinity, we have 

^<«'*<->i = < 45 ' 
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Figure 2: The partial energy absorption of the scalar wave with the change of 7, I 



and therefore we have the jost functions 



(±) 



r 2 VAB 
2ibjr 2 H 



(46) 



In numerical calculation, we would like to express gi in Eq. (|32[) and Eq. (|33[) in terms of the jost functions, and the 
absorption cross section therefore reads 



'tbs 



and 



abs 



TT 21+1 

^\ft ] \ 2 ' 



7T ^ 21+1 



(47) 



(48) 



We now report our numerical result for the massless scalar perturbation in the CFM braneworld black hole back- 
ground. Inserting the metric components into Eq. (|10p . the effective potential has the form 



V s ,i(r) 



1 



1(1 + 1) + 2M 



Af (7 - 3)(r 2 - 6Mr + 6M 2 ) 



r 3 (2r-3Af) 5 



(49) 



Taking M = 1, the bahavior of the potential is shown in Fig.l. We see that there is a potential barrier outside 
the horizon and for the fixed angular index I this barrier increases for more negative 7. Fixing 7, we find that the 
potential barrier increases with the increase of the angular index. When 7 becomes small, the potential has a negative 
well after the positive barrier. The negative well appears earlier for very negative 7. This was also observed in [To! ]. 
Comparing with the height of the positive barrier, the absolute value of the negative peak is very small for chosen 7. 
Thus the potential barrier will dominate in the determination of the wave dynamics, which qualitatively leads to the 
general quasinormal ringing as disclosed in [ic| . 
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Figure 3: The total energy absorption of the scalar wave. 



The spectrum of absorption cross sections for both Schwarzschild and the brane- world case with 7 = — 5 are shown 
in Fig. 2a for different I. Fig. 2b displays the partial absorption cross section for chosen I and different 7. It is noticed 
that the absorption cross section of the brane-world black hole can be either smaller or larger than that in the 
Schwarzschild black hole provided that 7 is smaller or bigger than 3. This result is consistent with the behavior of 
the potential barriers as shown in Fig.l. The higher or lower barrier of the potential could decrease or increase the 
absorption of the scalar field around the brane-world black hole. For the low u region, we observed that the absorption 
is higher for smaller values of 7. This could be the imprint of the negative well in the potential, since different from 
the barrier, the potential well enhances the absorption rate. But the overall behavior of the energy absorption shows 
that the potential barrier dominates the physics. 

In Fig. 3, we show the spectrum of the total absorption cross section, we see that with the decrease of 7, there are 
fewer oscillations. For smaller 7, in the low energy band the total absorptions are suppressed. This is because of the 
increase in the potential barrier which plays the dominate role for smaller 7 leading to the decrease of the absorption 
of the incoming wave. The full penetration occurs when the incoming wave with very high energy, which explains 
why the total absorption approaches a constant value at high frequency. 
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III. GRAVITATIONAL PERTURBATION AND ITS ENERGY ABSORPTION CROSS SECTION 

Now we start the formulation of gravitational perturbation, of which the detailed deduction requires special tech- 
niques provided in [25l . [29l | . We will concentrate our attention here on the brane-localized axial perturbation in the 
background of the braneworld black hole. 



A. Axial Gravitational Perturbation Equations 

Using the formalism introduced in [2^|, the axially perturbed metric Eq. (Q]) is 

ds 2 = e 2u {dtf - e 2lp (d V - wdt - q 2 dx 2 - q 3 dx 3 ) 2 - e 2 " 2 {dx 2 ) 2 - e 2 " 3 {dx 3 ) 2 . (50) 

Here 

x 2 =r, x 3 = 6, (51) 

w and qt (i — 2, 3) are quantities that embody the perturbation. Due to the lack of knowledge of the 5D bulk metric, 
we restrict to the brane localized perturbation governed by the perturbation equation 

SR^ = 0. (52) 

This simplification is supported by the analysis of gravitational shortcuts [30| which shows that gravitational fields do 
not travel deep into the bulk. It can be justified at least in a regime where the perturbation energy does not exceed 
the threshold of the Kaluza-Klein massive modes. By defining 



one finds equations 
from 8R\2 = 0, and 
from SRi3 = 0. We take 



Qij — qi,j — Qj,i , Qoi = —Qi,0 ! (53) 
( e 3^-M2-M3Q 23 ) )3 = _ e 3^+M3-WQ 02)0 (54) 

(e 3^-M 2 -„3Q 23))2 = e 3^+„3g 03)0 (55) 



e 3^+„- M 2-wQ 23 = r 2 VAB sin 3 9Q 23 = Q (56) 
in Eq. (|54[) and Eq. (|55|) . and then cancel w 1 to arrive at an equation of Q: 



2 d Q , 4 rrs 9 f v AB dQ \ , , . 3 „ d ( 1 dQ 



~ r + r v AB— ^--^ + Asm 3 9— — ^--2. = . (57) 

dt 2 dr \ r 2 dr J d9 \sm 3 9 09 J y ' 

Analogously, this equation plays the same role as the Klein-Gorden equation in the scalar perturbation. However the 
equation here cannot be separated in terms of Legendre polynomials as in the Klein-Gorden equation, the Gegenbauer 
polynomials are employed to carry out the separation [25j. Q is now seperated as 

oo 

Q(t,r,8)=Y / rZi(r)C- + 3 2 /2 (cos9)e-^ t , (58) 



1=2 



where the Gegenbauer polynomials C ;+ 3 2 ' /2 (0) satisfy the equation 



,2 VQ d {„--2vn dC "( COs( 



d§ l sin ~ — d9 — ' +1(1 + 2v)Ci{cob9) = 0. (59) 



The radial Schrodinger-like equation can be derived in the form, 



,hZl +[u; 2 -V g Ar)}Zi = 0. (60) 



2 



is the tortoise coordinate defined by dr* = drj\J AB, and the effective potential is 



M,) = A "- 1 »|' + 3 » + ^-^^*. (61) 
r r Ir 
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B. Boundary Conditions 

We set the boundary conditions to be those of plane wave scattering. 



1. Boundary Condition at Infinity 

It is expected that the quantity Q = r 2 V AB sin 3 9Q23 has the similar behavior of the scattered spherical wave since 
a straightforward examination shows that Q(t, r, 9) = e l (^ rcose -^ t ) [ s exactly a solution of Eq. (|5"7| in the region far 
from the black hole. Therefore it is appropriate to write 



Q(t,r,6) 



„iujr cost) 1 



(62) 



For the infinite r, Q behaves as Eq. (|62[) , which can also be expanded in a combination of Gagenbauer polynomials 
as the decomposition of Q in Eq. (f55j) . The expansion of the plane wave part is [3lj 



jiuir cos 1 



1 00 

= -o E ( l + i/2)^M^Mcr + 3 2 /2 (cos0). 



(63) 



Note that in practice we calculate only the partial absorption cross section of low angular indices, thus when ur is a 
large number, the asymptotic approximation for jz(wr) reads 



ji{x) ~ - sin(x ~ ln/2) = -±- [e" M - (- I 
x 2x 



l Ax~\ 



(64) 



for x ^> 1 and x 3> I, where x — wr. There exists a certain Iq and before this Iq the summation terms of Eq. 
may be replaced by the asymptotic form, 



Aujr cos 9 



IJ2(1 + l/2)(-l) V) [e— - (-1)^] C- 3 / 2 (cos( 



2=-2 



1 00 



(65) 



The larger lot is, the larger Iq will be. 

We also need the expansion of the term representing the scattered wave. Formally we may write 



f{6) c ikr 



E »MC;+2 /2 ( COS 



(66) 



(=-2 



Similar to the plane wave, we expect that before some integer Vq, the summation coefficients in Eq. (166j) are asymp- 
totically 



Xi(r)~ \{l + l/2){ur)f?j»\ 



(67) 



at some very large wr, where are complex numbers depending solely on index Z. Thus we can write the expansion 
Eq. UBBJ as 



/(<?) 



1 



e<W r = - £ (Z + l/2)M/f e-C- 3 2 /2 (cos ( 



l=-2 



+ E x^w^r+f^ose). 
i=i' +i 

Combining Eq. ((65)) and Eq. (|68j) . we have the applicable form of the asymptotic behavior of Q(t, r, 9): 

„-iwt miu (lo,l' ) 



(68) 



6 



(1 + l/2)(-l)'(wr) [e-^ r - (-l)'Sf e^ r ] C^ 2 {cos0) + 



(69) 
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Here the omitted terms are those over the upper limit for applying the asymptotic form, which are not applicable to 
our numerical calculation. Sf = 1 — jf is a complex number being the ratio of incoming and outgoing wave. The 
summation starts from index 2 for the modes from / = — 2 to / = 1 are not physical in gravitational perturbation. 
Since min(7o, Iq) certainly exceeds our need for numerical calculation, we may write for simplicity 

-iut 00 

Q~—J2(l + l/2)(-l)V) [e"^ r * " (-l)^V^*] C- 3 2 /2 (cos0), (70) 

1=2 

which is valid for a limited number of angular indices. As done in the scalar case, in the exponentials we have replaced 
r into r*, since in the remote region the spacetime is flat and we actually have dV* ~ dr. Therefore the boundary 
condition for Z\ (r) at infinity reads 

Z,(r)~i(l + l/2)(-l) , w[e- fa ' r '-(-l) , 5fe fa ' r *] r -> oc. (71) 



2. Boundary Condition at the Horizon 
We now study the boundary condition in the vicinity of the event horizon rjj. Expanding Zi(r) near the horizon [27| 

dZ L ^Pjp(j L _ 1 Y-\ *Zt„P«tp-l) fjL-xr', (73) 
dr r ff / dr 2 - r z H \r H J 

when r — ► rjy. = ci t o with the dimension of length -1 . Considering Eq. (|18[) . in the vicinity of the event horizon Eq. 
(|73|) becomes 

2r a u A / „ \ p 



we have 



r 



d?t \r ff / dr 2 r H \r H J dr r H \r H J 

and 

[co 2 -V eff (r)]Zi ~ PlU > 2 (j- - lj P . (75) 

Inserting Eq. (JHJ) and Eq. ([75]) into Eq. ((501), we obtain 

A 1 B 1 [p(p-l) + p}+r 2 H uj 2 = 0, (76) 

where 



p - ±iu>r„/ y /A 1 B 1 . (77) 

Similar to the scalar case, we choose the minus sign in the p expression since there is only ingoing wave near the event 
horizon and nothing can escape from the black hole. The boundary condition at the horizon then reads 

Zt(r) ~ Pl (r/r H -iy^ r »/^^, r^r H + 0. (78) 

In accordance the behavior of Q(i, r, 9) at the horizon is 

oo 

Q(t,r,6) ~ e- iut r H (r/r H ~ lP^/^™^^ 3 / 2 ^) , r -> r ff + 0. (79) 

i=2 



C. Application of the Newman-Penrose Formalism 



The study of the gravitational perturbation and the gravitational energy absorption spectrum involves considerable 
algebraic complexity. There is a well developed approach to this problem provided by Teukolsky [26|, where 
Newman-Penrose formalism was employed. In this section we will employ the formalism to calculate the energy flux 
coming from infinity and the energy falling into the black hole. 
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1. Formalism in the Kinnersley Tetrad 

This tetrad serves to calculate the incoming energy flux of the gravitational wave at infinity(32|. It is set to be 

i» = i(i,V3s,o,o), (80) 

n» = 1(1, --/AS, 0,0), (81) 

1 (O,O,l,icsc0), (82) 



V2r 

of which the components queue in the order (t, r, 9, ip). The corresponding spin coefficients are found to be 

K = a = X = v = e = TT = T = 0, (83) 
1 [B n cot 9 y/AB A r [B 

P= rVA' a = - p= 2V2>' M= ^' 7= fVl" (84) 

Here, taking A = recovers the Schwarzschild case. 

Then we need to evaluate two contractions of the perturbation of the Riemann tensor with the Kinnersley tetrad, 
which play key roles in calculating the energy flux at infinity and the energy falling down the hole, and which hereafter 
is denoted by and \ defined as 

¥^ = -5R^ Xa l»m»l x m° \ (85) 

and 

#W = -6R^x<jn f *m"*n x m a *. (86) 
To perform the contraction a new orthonormal tetrad was introduced 

e^ 0) = (e- v , 0,0,0), (87) 



ft, = (O.e-^0,0), (88) 



e (i) 

e^j = (0,0,e^ 2 ,0), (89) 

ef 3) = (0,0, 0,e-^), (90) 

where components are ordered by (x , x , x 2 , x 3 ) = (i, <p, r, #). The components of Kinnersley tetrad against this 
tetrad basis are therefore 

= (e-",0,e- v ,0), (91) 

n W = i(e w ,0,-e«',0) ) (92) 

m<"> = -L(0,i,0,l). (93) 

Here the indices in the bracket indicate the components against the basis of Eq. (|57|) through Eq. ([9"U]). and they run 
through to 3. Now we do the contraction of Eq. (|85)) and find 

*o 1} = - ie" 2ly ((5i?030i + SR2321 + SR2301 + SRo32i) 

~2 e 2l, (^0303 + 2SR0323 + SR2323 ~ SRqiqi ~ 2SR i2i — SR2121), (94) 



as well as 



= ~ ^e 2ly (-(5i?0301 - SR 23 21 + SR2301 + SR0321) 

~ o e2iy (^^0303 — 2(Si?Q323 + ^2323 — ^-RoiOl + 25i?Q121 _ <^212l)- (95) 
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Here, in both expressions, the second bracket would vanish in purely axial perturbation, whose demonstration was 
raised in §31 of (25j and applicable to all static, spherically symmetric metric. For convenience, we denote the non- 
vanishing part, which contains the first bracket, of Eq. (|94| and Eq. (|95f to be i Im^g 1 ^ and i Im^4 , and the 
vanishing part, Re*& and Re^ . The Riemann components therein can be simplified as (§13 of [25]), 



<5-Ro30i 

<5^2301 
<^0321 



1 



,4>-2lL 2 -V3l 



lp-2fJ,2-(J, 3 



tp—V—IJ,2-IJ, 3 



7TM2,r J <?23 + -zQ23,r 

r I 1 



1 



1p—U—IJ,2—ti3 



<?03,2 — <3o2,3 

20,3 + Q20 COt 9 — Q23,0 + Q30 



- Q02 cot 

1 



V 



Inserting these equations into Eq. and Eq. 



we have 



sin 9 



lff (i) _ O23 , 1 



2<^23,r + 2 6 



Qm,r 



\l2,r Q 



1 



2i/ ) 2/i2 /O 
e V03,0, 



and 



4 e -2i/+2/i 2 



smt 



j ,f,(l) ^23 J-„ 

Im^ = 2 ( 5 23 ' r + 2 6 



' ^ +Al2 Q 23 ,0 + i(M2,r+^)Q 



23 



— e"^ 2 



/03,r 



f,r + M2,r ) <?03 



l e -2 I .+2 M2( 



;o3,o- 



With the aid of Eq. ([54]) through Eq. (|57|) . the two contractions can be simplified to be 



2iu Im* 



(i) _ 
— 



E 

;=2 



2v/AB A'B + AS' 



2V AB 



and 



8iw Im* 



(i) 



E 



2VAB A'B + AB' 



2VAB 



— 2iuj A.r\Zi 



2iu; K w Zi + 



iu>(AB' - A'B) 



2VAB 



iu}(A'B - AB') 



V, 



9,1 



C; + 3 2 /2 ( COS( 



rA 2 sin 2 1 



2VAB 



+ V9.1 



C i+ 3 2 /2 (cos( 
r sin 2 9 



(96) 
(97) 

(98) 

(99) 



(100) 



(101) 



(102) 



(103) 



where V 3t i is the radial effective potential defined by Eq. (|6~Tj) and \±) = d/dr* ±iui. It is easy to check that equating 
A(r) and B(r) would recover the result in the Schwarzschild case. 



2. Formalism in the Hawking-Hartle Tetrad 



To calculate the gravitational wave energy flux on the horizon, the Kinnersley tetrad is no longer applicable because 
of their singular behavior at the horizon. We use the tetrad introduced by Hartle and Hawking instead[33|], which 
represents a physical observer crossing the event horizon. It eliminates the singular behavior by imposing on the 
Kinnersley tetrad a rotation of the third class ([2||, §8). In our problem the rotation parameter is A = 2 /A. The 
tetrad is normalized so that on the basis of coordinate (v = t + r*,r,6, tp), {l HH ) v — 1. This results in the Hawking- 
Hartle basis as follows: 



(i HH r 

(n HH r 



1/2(1, VAB, 0,0), 

(i/a,-/b7I,o,o), 



(104) 
(105) 



and is unchanged. The components are arranged in the order of (t, r, 9, (p). Changing to the basis of (d v ,d r , dg, d v ), 
we have 



(i HH Y 

{n HH Y 



(i,Vab/2,o,o), 
(o,-v^b7I,o,o), 



(106) 
(107) 
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again with unchanged. Thus we obtain the well-behaved tetrad on the horizon, which is equivalent to a physical 
observer advancing along the direction l HH . Hereafter we perform calculations in this Hawking-Hartle tetrad. We 
can derive new spin coefficients using their transformation rules under the third class of tetrad rotation ([25j], §8), the 
one of importance is 



HH 




where D is the differential operator corresponding to the tangent vector We define 



co 



HH 



1 

4 



<B dA 
A dr 



(108) 



(109) 



2eo is the surface gravity of the black hole ([25|], §8). We also need to know that n vanishes globally, ensuring the 
integral curves of l HH to be null geodesies, and that p HH vanishes at the horizon, since 



a-v 



IAB 
2r 



0, (r^r H ). 



(110) 



Moreover, we present here two relations needed for calculating the energy falling into the hole. They are obtained 
with the perturbation effect, where the Hawking-Hartle tetrad is also perturbed so that l HH is the generator of a 
congruence of null geodesies, and n always stays zero. We start from two Ricci identities under any Newmann-Penrose 
basis, which are 



Da — Sk = cr(3e — e* + p + p*) + k(tt* — r — 3/3 — a*) + 4"0) 
Dp-5*n = (p 2 + \a\ 2 ) + p(e + e*) - h*t - n(3a + f3* - tt) + $ o 



(111) 
(112) 



f [25j . §8). Here \l/o and $00 are projections of the Weyl tensor and the Ricci tensor upon the tetrad. Linearizing these 

s(l) 



two and with the consequence from 5R^ = that %y = 0, we find, under Hawking-Hartle tetrad, 

D HH p HH(l) = 2eQp HH(l)_ 



(113) 
(114) 



D. Gravitational Wave Energy Absorption 

1. Incoming Energy Flux at Infinity 



The general expression for the incoming energy flux is [25l 1281 13 



d 2 E^ 



dtdn Mttlu 2 r 



' lim r 2 !^ 1 ^ 2 , 



(115) 



where is defined in Eq. ([94|). Considering the axial perturbation, we can replace $q by Im^Q 1 ''. Taking r — > 00 
in Eq. (|102p and taking into account only the incoming wave, we have 



OO 



1=2 



rsin 9 



(r — > 00). 



With the boundary condition at infinity Eq. (|71|) . we have 



A ( _)Z; (mc) — ► (-iy +L iuj(2l + l)e~ wr * /6, (r -> 00). 



(116) 



(117) 



Inserting Eq. (|116p . Eq. pi7p into Eq. (| 1 1 5[) and integrating the flux over all directions, we obtain the total incoming 
energy per unit time at infinity 



(£(inc) 



in 



dtdQ. 



,2 00 



2304tt 



E 

1=2 



in 



Cl+2 2 ( COSt 

sin 2 6» 



to 

192' 



(118) 



15 



The last equality was obtained by using the orthogonality of the Gegenbauer polynomials [3l| 

1 ^-3/2^-3/2 



and the summation 



dx = (119) 

(1-x 2 ) 2 (2/ + l)(i + 2)(/ + l)i(I-l)' 1 j 



V 21 + 1 - 1 

^(Z + 2)(Z + l)/(/-l) _ 3- (1 } 



2. Energy Falling into the Hole 

Calculating the gravitational wave energy absorption into the black hole is more complicated than that in the scalar 
case. We will follow Teukolsky's treatment [28| to deal with this problem. 

With Eq. (|114p . following the procedure presented in [2^|, §98, we have the change in the area of the event horizon 
in the process of perturbation 

i!^_ii k ^(i)| 2j (121) 



dvdQ, e ™ 

where S denotes the area of event horizon, and the spin coefficient eo is presented in Eq. p09p . The first law of black 
hole thermodynamics relates the change of the horizon area to the change of the black hole internal energy 

dE = — dS, (122) 
8tt k j 

where R = 2eo is the surface gravity. Thus 

d 2 S 4tt d 2 E^ 



dvdQ eo dvdQ 

Combining Eq. (fT22| and Eq. (JT23)), we have 



(123) 



The spin coefficient a^J 1 ^ is evaluated in Eq. (| 1 13|) . where 

D?J* = d v =d t ^ -iw, (125) 

and consequently we get 

a HH(l) = _*0 U (126) 

H iuj + 2e 

Referring to the transformation of Weyl scalars under the rotation of the third class ((25|, §8), we have 

y" H{1) =A 2 'f { 1) /4. (127) 
Combining Eq. ([T2T|) . Eq. (fT26|) and Eq. (fT24"]) , we arrive at 

d 2 E (abs) ^ \A 2 ¥ 1] \ 2 rH 



_ [ (128) 

dvdn 16tt4w 2 + A x B x r- 2 ' K ' 

Now we need to evaluate = —SC fl „ p \l tl m l/ l p m x 1 

= - ie- 2l, (6C 30i + SC2321 + 6C2301 + SC 032 i) 

--e~ 2ly (£Co303 + 25C 323 + SC 2 323 ~ SC 10l - 28C 121 - SC 2 12l)- (129) 
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The Weyl tensors are defined as 



C^ U p\ — Rfj, V px — —(g^pRxv — g^xRpis + gu P R\ P — gu\R PP ) + -^R(g PP g\v — g P \gpv)- 



(130) 



Three facts can help greatly simplify the evaluation, which are, (i) only the diagonal components of the Ricci tensor 
is non-zero for a metric like Eq. |T]), (ii) according to the perturbation equation, SR^ = 0, and (iii) only 5g\^ 7^ 0, 

(p = 0, 2, 3). The evaluation shows that the components of the perturbed Weyl tensors in vI/q 1 ^ all coincide with those 
of the Riemann's, except for <5Co3oi and SC2321, which are 



SC. 



0301 



SR, 



0301 



RooSg?. 



31 



<5C2321 — <5i?2321 + 2^22^331 



Therefore 



(i) 



(i) 



-2i/ 



■(Rao + R22)Sg 3 



1 



-RgsaSg 



31- 



-2// 



-R(goo + 522)^331 ■ 



(131) 
(132) 

(133) 



This is further simplified by the fact that we carry out the contraction under the orthonormal basis, Eq. (|87|) 
through Eq. ([^0]) . where for a spherical metric as Eq. (Q]) describes, there holds goo = ~.922 = 1, implying that 
(goo + 522)^531 = 0. And also, we have the extended result from the calculation of the Schwarzschild model 



Roo — — 



BA" A'B' 



BA> ' — (A') 2 - 



R22 = 



2A AA rA 
BA" A'B 1 B' B , 



2A 



AA 



which gives that 



R 



00 



R22 = ^(^-^)~0(r/r H -l), (r^r H + Q). 



(134) 



(135) 



That is, the term Rqq + R22 is of the first order small quantity as r tends to the horizon. On the other hand, it is 
easy to check that near the horizon, we have (73,2 ~ Q(r/rH — 1) _1 ~ {r/r^ — l^-K^/v^TSij-l^ so i ag ^ $g 31 ^ g 3 ~ 
(r/rii — \)^ l ^/^ AlBl \ and hence we have (Roo + i?22)£<?3i — > 0, as r — > r# + 0. Therefore in our context 



= i lm¥ Q 1] . 



Hence 



d 2 E ( - a 



hs) 



r\ \A 2 \nx^ 



(1)|2 



dvdtl 16-7T Acu 2 + A 1 Bir H 2 ' 
In the vicinity of the event horizon, r — > m + 0, applying the boundary condition Eq. (1781) . we get 

cr + 3 2 /2 (cos0) 



(136) 



(137) 



1=2 



2iuirH sin 9 



00 

[(/MBlr^ - 2iw) pi(r/r H - i)-^r H /VMBl 



1=2 



C- 3 2 /2 (cos0) 
rn sin 2 9 



(138) 



Inserting this expression into Eq. (| 1 3T|) and performing the integration over all directions, we have the total absorbed 
energy per unit time 



^2j^(abs) 

dvdQ 



dQ 



\_Pll 
16tt 



C- + T(cos9) 
sin 2 9 



Employing the orthogonality of the Gagenbauer polynomials, Eq. (138) becomes 



N 2 q-2)! 

i^2l+l(l + 2)V 



(139) 



(140) 
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3. The Absorption Cross Section 

We define the absorption cross section as the ratio of the total absorbed energy to the total incoming energy at 
infinity. Averaging the total incident energy over the event horizon, we have the incident energy flux in the analogous 
sense, that 



The total absorption cross section reads 



3 4nrl 1536^ 



i=2 

The partial absorption cross section has the form 



_ $W 1728^ ^ h| 2 {l-n (U9] 
9 j(™c) w 2 Z^2l + l(l + 2)V [ ' 



ahs 

a 



^ 2/ + l(/ + 2)i- (W3) 

Now we proceed to evaluate pi, which is necessary to obtain the 

By means of the Wronskian of the radial perturbation equation Eq. ([50|) . for any two solutions and Z?\ their 
Wronskian is 

W[Z?\ZW} = ^=, (144) 

where c| 12 ^ is a complex number depending on the explicit form of Zj 1 ^ and ZY , Applying this to the asymptotic 
solution both at infinity, as is shown in Eq. (|7ip . and at the horizon, as is shown in Eq. (178|) . we have, for the region 
at infinity, 



iuj 3 (2l + l) 2 , , ol9 , iu 3 (2l + l) 2 „ , s 



and for the region approaching the horizon, 



Equating the above two equations gives 



where 



W[Zt,Z l ] = -^=\p l \\ (146) 
V j4£> 



|p , |2 = ^Z rf , 



rf = i - isfi 2 , (148) 

Inserting this |p;| 2 into the expression of partial absorption cross sections, Eq. (|143[) . we get 

a^ = 12^(2? + l)|^||rf, (149) 

which is analogous to Eq. (|38[) in the scalar case. In the high energy limit, the wave can easily penetrate the potential 
barrier and be absorbed by the black hole, Tf — ► 1. Thus summing up both sides of Eq. (|149j) and applying the 
summation formula Eq. (|120[) . we get 

a abs — >47rr 2 f) ( w _>oo). (150) 
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E. Conservation of Energy 

Here we will justify the formalism for the energy absorption rate in §III-D-2. In the far away region, we can write 
Eq(70) in the form 



Q ~ e ~ wt J2 r (Xe^ wr * + K t e lujr *) C~ 3 / 2 (cos 0) , r -> oo. (151) 



-3/2 

1=2 

Near the horizon we can follow Eq. (|79p and express 



THir/TH-ir^WXttYtWi+i i( cos 0)> r-r* + 0, (152) 



Z=2 

where we have replaced p; with 7; to avoid ambiguity. Our goal is to show that the difference between the incoming 
energy and the reflected energy by the potential barrier at infinity is all absorbed by the black hole 

dE (abs) ^ d£(mc) ^(re/) 

c?t fit rft 

For the incoming energy per unit time, we adopt Eq. 1)1 15p . and get 



£§£h„ ,154, 



dt ^ 16tt 

1=2 



where we have employed the orthogonality of the Gegengauer polynomials and we let 

> 

dn = Ei (155) 



C /+2 /2 ( COSr ' 



sin 2 6 

for simplicity. The reflected energy can be obtained similarly. Referring to ref. [25|, §98, the out flowing energy at 
infinity is 

Inserting Eq. (|103|) into it and using Eq. (|15ip . we have 

~ |7?, 1 2 

< 157 » 

/=2 

The absorbed energy per unit time can be got by referring to the calculation in §III-D-2 and using Eq. (|152p 



di ^ 16tt v ; 

1=2 

To complete the demonstration, we need to establish the relation among X;, IZi and T\. We do this by means of 
Wronskians of the radial equation Eq. (|60|) . Referring to the last section where we have calculated the Wronskian, 
and by virtue of Eq. (|151[) and Eq. (|152[) . we have the Wronskian at infinity that 

W\Z* l ,Z l \ = 2iu (-|X| 2 + |7^| 2 ) {AB)- 1 ' 2 , (159) 

and also the Wronskian at the horizon that 

W[Zf,Zi] = -21^ {AB)- 1 / 2 . (160) 

Equating the above two, we have the relation among X;, IZi and X/ that 

\T l \ 2 = \I l \ 2 -\K l \ 2 . (161) 

Now we can relate Eq. (|154[) , Eq. ()157|) and Eq. ()158|) by using Eq. (jl53[) , thus completing the demonstration. This 
demonstration shows that the application of Hawking-Hartle tetrad in our above discussions is well grounded on the 
law of energy conservation. 
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Figure 4: The effective potentials of the gravitational perturbation for different 7 and I. 

F. Numerical Study for Gravitational Perturbation 

In solving the radial equation Eq. (|60[) . we choose a particular solution 4>i(r) to be normalized as 

4> l {r)=pJ 1 Z l {r). (162) 
In the remote region, it consists of the outgoing and ingoing waves, which is 

fa (r) = n\- ] M l(+) + n\ +) Af K - h (163) 

where 

00 

M(±)(r) = ^- E"£ } (t) B ' (164) 
e called jost i 

their Wronskian, by applying Eq. (|144[) we arrive at 



71=0 



and the combination coefficients n ; are called jost functions. To obtain the jost functions numerically, we calculate 



The jost functions read 



WW K+) ,M(-)} = -S=. (165) 



», {±) = T^W[M {±) }. (166) 
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Figure 5: Partial energy absorption of the gravitational wave for different 7,/. 



Comparing Eq. (|163p with Eq. (fTTj) . we have 



Pi 



s? 



v (2[+l> 
[ ' - (-) ' 



12 n 



1 

(+) 



(-1) 



1+1 'n 



,(-)' 



(167) 
(168) 



They are key quantities for determining the absorption cross section. 
In practice, we insert Eq. (|167|) into Eq. Ijl42|) and Eq. (I143|) . and have 



'9,1 



= Ylur 



(l-2)\ 21 + 1 
(l + 2)\\n\-Y 



and 



19 2 2^ + 1 

12 ^2. (/ + 2)!|n (- )l2 - 



;=2 



(169) 



(170) 



Now we can do the numerical calculation. Substituting the CFM black hole metric, for the gravitational perturba- 
tion, the effective potential has the form 



V g ,i(r) 



1 



1(1 + 1) 



6M 



M(7 - 3)(5r 2 - 20Mr + 18M 2 ) 



r 3 (2r- 3M) 2 



(171) 



Choosing M — 1, we display the behavior of the potential for the axial perturbation in Fig. 4. The effective potential 
is not positive definite. For enough negative value of 7, a negative peak will show up out of the horizon. The negative 
potential well becomes deeper when 7 becomes more negative for fixed I. For chosen 7, the negative potential well 
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Figure 6: Total energy absorption of the gravitational wave. 

appears for small /. Compared with the scalar case, in the gravitational perturbation, the negative wells appear before 
the potential barriers. It was argued that negative potential may result in the amplification of the perturbation out 
of the black hole and cause the spacetime to be unstable [35[. However in the QNM study of the CFM braneworld 
black hole [l(|, it was shown that even for very negative value of 7, the perturbative dynamics out of the black 
hole is always stable. The negative potential well loses the competition with the positive potential barrier and the 
perturbative dynamics is still dominated by the positive potential barrier. 

In the study of the absorption, the potential well tends to enhance the absorption while the barrier tends to decrease 
it. The result on the partial absorption cross section is shown in Fig. 5. We find that as the case observed in the QNM 
study (lo| , the potential barrier wins the competition with the well and dominates the contribution to the absorption 
rate. Fixing I, we observe that for more negative 7, less gravitational wave energy is absorbed by the CFM braneworld 
black hole. And for chosen 7, the absorption rate decreases with the increase of I due to higher potential barrier. In 
the low lo region, the absorption is enhanced for smaller values of 7, which is the effect of the negative potential well. 

In the high energy limit, the partial absorption cross sections of the gravitational wave flatten out which is different 
from the case in the scalar wave. This is due to the difference in Eq. (|38|) and Eq. (|149[) . where in the scalar case the 
partial section is proportional to lu~ 2 in the high frequency regime, while in the gravitational wave case, the partial 
section is proportional to oj° in the high frequency regime. Deep reason causing this difference lies in the different 
expansions of the plane wave. In the scalar case, the plane wave is expanded in terms of Legendre polynomials into 
spherical waves with a factor (ojr)~ 1 e ± ' lujr (Eq. [12]), while in the gravitational case, it is expanded in Gegenbauer 
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polynomials into the wave with a factor (uir)e ± ' lujr (Eq. [65]). 

The result on the total absorption cross section of the gravitational wave is shown in Fig. 6. The total absorption 
cross section, which is the physically observable quantity behaves in a similar manner to that in the scalar case. This 
shows that the difference in the gravitational partial absorption cross section from that of the scalar case is more 
mathematical than physical. 



IV. CONCLUSION AND DISCUSSION 



In this paper, we have studied the energy absorption problem of the scalar wave as well as the axial gravitational 
wave in the background of the CFM brane-world black hole. The CFM black hole is spherical and has only one event 
horizon at rn = 2A1. When the parameter 7 = 3, the CFM braneworld black hole returns to the Schwarzschild 
solution. In comparison with Schwarzschild black hole, the CFM braneworld black hole will be either hotter or colder 
depending upon whether 7 < 3 or 7 > 3. 

We have calculated the scalar perturbation and the axial gravitational perturbation around the brane black hole 
and evaluated the energy absorptions of the scalar and gravitational waves. Comparing with the 4D Schwarzschild 
black hole when 7 = 3, we observed that the energy absorption of the braneworld black hole decreases with the 
decrease of 7 starting from 3. While when 7 > 3 we found that the energy absorption enhanced compared with that 
of the Schwarzschild black hole. We restricted to the case 7 < 11/3. This result holds the same for both the scalar 
wave and axial gravitational wave outside the CFM hole. In both perturbations, the negative potential well appeared, 
however the positive potential barrier still dominated in determining the absorption rate, which has the same effect 
as observed in the QNM study. The result on the absorption spectrum implies that for 7 < 3, the black hole emission 
will be enhanced with the decrease of 7, while for 7 > 3, the emission of the black hole will be suppressed, which is 
consistent with the behavior of the Hawking temperature of the braneworld black hole. We conclude that the energy 
absorption for the scalar and axial gravitational wave gives signatures of the bulk effects in the brane-world black 
hole, which differs from that of the 4D Schwarzschild black hole. We expect that these signatures can be observed in 
the future experiments, which could help us learn the properties of the extra dimensions. 

In the study of the axial gravitation perturbation, the deduction of = helps us a lot in doing the 

calculation, however the simplification is obviously not hold for polar perturbation. Thus it is of interest to generalize 
our discussion in the future to study on the polar gravitation wave. 

It needs to be emphasized that although several spherically symmetric and static brane black hole solutions with 
contributions from the bulk gravity have been found, none of these are obtained as exact solutions of the full five- 
dimensional bulk field equations [15[. The propagation of gravity into the bulk does not permit the treatment of the 
brane gravitational field equations as a closed system. This fact limits the consideration of propagating modes of 
fields off the brane. Although the propagating modes in the bulk are hard to be obtained at the present moment, the 
modes on the brane are still interesting since they are the most phenomenologically interesting effects which can be 
detected during experiments. Furthermore, it was argued that the emission of particle modes on the brane is dominant 
compared to that off the brane Brane localized modes have been investigated to disclose the information of extra- 
dimensions in different attempts [id [Til . 0, [l?], 0, Oil ■ Despite that we cannot obtain the propagating modes in the 
bulk due to the lack of complete bulk solutions owing to the conceptually complicated gravitational field equations, 
our study of brane modes is still well-motivated and interesting. Employing the Newman-Penrose formalism, we have 
provided a test of the energy absorption spectrum of a new braneworld black hole solution. 
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